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ESTIMATION  OF  BLADE  AIRLOADS 
FROM  ROTOR  BLADE  BENDING  MOMENTS 

William  G.  Bousman 
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U.S.  Army  Aviation  Research  and  Technology  Activity 
Ames  Research  Center 
Moffett  Field,  California  94035 


1.  ABSTRACT 

A  method  is  developed  to  estimate  the  blade  normal  airloads  by  using  measured 
flap  bending  moments;  that  is,  the  rotor  blade  is  used  as  a  force  balance.  The  blade’s  ro¬ 
tating,  in  vacuum  modes  are  calculated  and  the  airloads  are  then  expressed  as  an  algebraic 
sum  of  the  mode  shapes,  modal  amplitudes,  mass  distribution,  and  frequency  properties. 
The  modal  amplitudes  are  identified  from  the  blade  bending  moments  using  the  Strain 
Pattern  Analysis  Method.  The  application  of  the  method  is  examined  using  simulated  flap 
bending  moment  data  that  have  been  calculated  from  measured  airloads  for  a  full-scale 
rotor  in  a  wind  tunnel.  The  estimated  airloads  are  compared  with  the  wind  tunnel  mea¬ 
surements.  The  effects  of  the  number  of  measurements,  the  number  of  modes,  and  errors 
in  the  measurements  and  the  blade  properties  are  examined,  and  the  method  is  shown  to 
be  robust. 


2.  INTRODUCTION 

Daughaday  and  Kline  [lj  described  an  approach  to  determining  generalized  forces 
from  measured  bending  moment  data  that  assumes  that,  at  each  harmonic  of  rotor  speed, 
the  blade  response  is  caused  by  the  blade  mode  nearest  in  frequency  to  the  harmonic. 
They  applied  this  single-mode  approach  to  flight  data  obtained  on  an  H-5  helicopter  with 
three  sets  of  rotor  blades  of  different  stiffnesses,  each  of  which  were  tested  at  three  rotor 
speeds.  They  calculated  the  generalized  force  for  each  harmonic  as  a  function  of  advance 
ratio  and,  although  there  was  considerable  scatter  in  the  flight  test  data,  they  were  able 
to  demonstrate  that  the  downwash  theories  that  w'ere  available  at  that  time  were  not  able 
to  predict  the  generalized  forces.  This  was  particularly  noticeable  at  advance  ratios  below 
0.20.  They  proposed  that,  to  the  extent  that  the  generalized  forces  are  independent  of 
the  rotor  blade  design,  the  generalized  forces  determined  from  flight  test  will  provide  a 
more  accurate  estimate  of  blade  loading  than  will  aerodynamic  theory  in  the  design  of  new 
rotors. 

The  possibility  of  using  generalized  forces  which  have  been  determined  from  flight 
measurements  on  a  helicopter  rotor  to  design  a  second  rotor  system  was  examined  in  more 
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detail  b>  Loewy  et  al.  [2j.  They  expanded  on  the  approach  of  Daughaday  and  Kline  to 
include  an  est  imat  ion  of  the  generalized  forces  using  the  three  nearest  blade  modes  instead 
of  a  single  mode.  They  calculated  the  generalized  forces  for  the  front  rotor  of  an  H  21 
helicopter  with  wooden  blades  and  applied  these  generalized  forces  to  an  H  21  helicopter 
with  metal  blades  and  to  a  YH-16  helicopter.  The  calculated  bending  moment  data  were 
then  compared  with  available  flight  test  measurements.  The  agreement  was  fairly  good 
for  the  rpetal- blade  H  21  rotor,  but  was  less  satisfactory  for  the  YH-16  rotor  Loewy  et 
al.  concluded  that  despite  the  differences  that  were  seen,  there  were  so  many  similarities 
in  the  generalized  forces  for  even  disparate  rotor  systems  and  flight  conditions,  that  the 
program  of  research  should  be  pursued.  They  also  concluded  that  the  use  of  three  adjacent 
modes  was  an  improvement  over  the  single  mode  approach. 

The  advent  of  flight  test  data  that  included  both  surface  pressure  and  bending  mo¬ 
ment  measurements  3]  encouraged  DuWaldt  and  Statler  to  re-examine  the  use  of  bending 
moment  data  to  estimate  the  generalized  forces  [4j.  They  treated  the  cosine  and  sine 
harmonics  separately,  and  in  this  way,  properly  accounted  for  phase  relationships  which 
had  not  been  done  in  t.he  previous  work.  They  compared  the  rotor-tip  deflection  that  was 
calculated  using  the  measured  airloads  with  that  obtained  using  the  measured  bending  mo 
merit.  data  for  one  flight  case  from  Reference  3.  They  concluded  that  the  comparison  was 
not  satisfactory,  primarily  because  the  methods  used  could  not  account  for  measurement 
errors 


Fscuber  and  Bousman  [5 j  have  recently  examined  the  structural  response  of  the 
CII-3;  rotor  to  the  aerodynamic  forces  1  fiat  were  measured  in  flight  and  wind  tun  tie! 
tests.  By  using  the  measured  airloads  they  have  been  able  to  eliminate  uncertainties 
in  the  aerodynamic*  model  and  in  this  way  study  the  structural  response  as  an  isolated 
problem.  One  result,  of  their  investigation  has  been  to  show  that  the  measurements  of  the* 
airloads  and  the  bending  moments  that  have  been  obtained  with  the  CH  34  rotor  are  quite 
accurate.  An  example  of  this  is  shown  in  Figure  1  where  the  3  to  10  harmonics  of  the  flap 
bonding  moment  measured  in  the  wind  tunnel  :6!  are  plotted  in  the  manner  of  Hooper  7 
using  a  three  dimensional  Cartesian  surface.  These  are  compared  to  the  moments  that 
have  been  calcuiat.-d  using  the  measured  a  reloads.  The  uncoupled  blade  flappmg  equation 
-  used  .or  the  ca'cu'ation.  To  obtain  the  good  agreement  that  is  shown  here  it  is  necessary 


t !  a» .  i  ' !  ■  ‘dude  si-  rfac  e  pressure  measurements  he  accurate,  2)  the  blade  flap  equal  ion  ol 
me.: ion  !•<:  ■. o r r *-c  t.  Ti  1  he  blade  mass  and  stiffness  properties  be  properly  comimfod.  4 }  tin 
solu’roit  alg-'riteni  .jsei  to  solve  1  be  differential  equations  be  sufficiently  a<  curate,  and 
•r>)  *  be  1  iade  bene! mq  moment  measurements  be  correct.  The  good  results  that  have  bet  n 
achieved  with  the  lor-e.j  response  cab  u'at  i- m  shown  in  Figure  1  raise  the  question  as  to 
vc1’-*!  cm  oi  ->t  o  w  i, •>•-.  ii>|e  to  work  ba<  k wards  from  the  bending  moment  measurements 
and  caV'dai*  1 m  spat  w  is»  airload  distribut  ions:  that  is,  use  the  rotor  blade  as  a  force 
bah. ;ic-’.  The  ?ior<  dynamic  loading  of  belie  of  er  rotors  is  still  not  well  uncic  rsiood  and 
1  t.e  a.uii,,  to  ohta.M  ;;  (too«{  wsr  up. ale  oi  ise  mrttia!  airload  c.  mg  b 

r  ’m»r  ton  vo i  id  provide  I  h<  <■:>. per  m.  ’  taust  a  sirnpie  and  inexpensive  technique  to 
m. '  a  oil  (to  or  .•  i  •  o'  v  i,  a  r  !<  '  a'  e  a  i  ’  \>  ot:  h  i  lv  es;  >ee uni  !y  valuaidc  ,!  t !  >•  •:.!!;•  >:  !* 
am;  r  >  inf  icm  ot  m  akmg  accurate  -.nri.o  ••  pressure  measurements  could  fa  avoided 


The  purpose  of  the  present  paper  is  to  demonstrate  such  a  technique  for  the  estimation  of 
the  normal  airloads  of  a  helicopter  blade. 


3.  FORMULATION 
Equation  of  Motion 

The  uncoupled  flapping  equation  of  Houbolt  and  Brooks  jS j  is  used  to  represent 
the  blade  dynamic  behavior 

(Elm")"-  (Tw')'  +  =F,  (]) 

where  El ,  the  blade  stiffness,  T,  the  tension,  and  m,  the  mass  are  all  functions  of  the  blade 
radial  coordinate.  The  blade  displacement,  rv,  and  the  normal  force,  Fz,  are  functions  of 
time  as  well  as  radius.  A  series  solution  for  the  blade  displacement  is  assumed  which 
separates  the  time  and  space  variables 

p 

™(M)  =  X^9"*^05^  +  ^nkasinknt)<(>n{r)  (2) 

n=  1 

where  qnkc  and  <7nka  are  the  amplitudes  of  the  kth  cosine  and  sine  harmonics  of  the  nth 
term  respectively,  <j>n{r)  is  the  nth  displacement  shape  function,  and  fl  is  the  rotational 
speed.  The  airload  distribution  is  separated  in  a  similar  fashion 

p 

Fz{r,t)  =  Yl^Fzkc^cosk^lt  +  Fzks{r)sinknt\  (3) 

n=  1 

where  Fzkc{r)  and  Fzks{r)  are  the  cosine  and  sine  harmonic  airloads  and  vary  with  the 
radius.  The  solution  is  then  expressed  for  each  cosine  and  sine  harmonic  separately 

p 

X;  qnkcKEItW  -  {T<f>'ny  -  mk2tl2<j>n I  -  Fzkc  (.1) 

fl"  1 

V 

X;  <7 nks\{ElK)"  (TKY  -  mk2 n2K\  -  Fzke  (5) 

n—  1 

As  has  been  shown  in  Reference  4,  if  the  series  chosen  to  solve  Eq.  (1)  is  assumed  to  be 
the  solution  of  the  homogeneous  form  of  that  equation,  that  is,  the  rotating  natural  modes 
of  the  blade  in  a  vacuum 


(6) 


WM)  =  7.  <t>n(r)sinunt 

n  —  1 

where  u>n  is  the  mode  natural  frequency;  then  substituting  into  the  homogeneous  form  of 
Eq.  (1)  gives 

Y^\{EI4>:r  -  (WJ  -  £  ra<t>nul  (7) 

n= 1  n—1 

and  eq.  (*l)  and  (5)  become 

p 

E zkc  -  ^  (w2  -  k 2ft1)m4>nqnkc  (8) 

n  -  1 

P 

F'zke  *  -  k2n2)m6nqnk,  (9) 

n  1 

where  the  order  has  been  reversed  to  show  the  dependency  of  the  airload  distribution  on 
the  modal  properties  and  modal  amplitudes. 

The  expressions  for  the  spanwise  harmonic  airload  given  by  Eq.  (8)  and  (9)  are 
useful  in  that  the  derivative  terms  of  Eq.  (4)  and  (5)  have  been  eliminated.  As  the 
number  of  modes  in  the  summation  in  Eq.  (8)  and  (9)  is  increased,  the  frequency  parameter 
(ou2  k“M2)  will  also  increase  for  n  greater  than  two  or  three.  At  the  same  time,  the  modal 
amplitudes,  qnkc  and  qnk«,  will  decrease.  For  the  summation  to  converge  it  is  necessary 
that  the  products  (ij2  -  k2M2)qnkc  and  -  k2M2)qnk»  become  small.  The  number  of 
modes  that  are  required  for  convergence  will  depend  on  how  quickly  this  product  term 
decreases. 

The  radial  variation  in  the  airload  will  be  determined  by  the  mass  distribution 
and  the  mode  shape.  The  frequency  parameter  (a;2  k2U2)  will  not  vary  with  blade 

radius,  nor  will  rt  vary  with  {light  condition  if  the  rotor  speed  is  constant.  The  only 
parameters  that  vary  with  fiight  condition  are  the  modal  amplitudes,  qnkc  and  qnk«<  which 
must  be  identified  from  measurements.  To  calculate  the  spanwise  airloads  then,  it  is  first 
necessary  to  obtain  the  mode*-  of  the  rotor  in  a  vat  uum  and  secondly  to  estimate  the 
modal  amplitudes  trom  blade  bending  moment,  data 

Modal  f’roper’te. s 

I  he  modal  properties  need  be  obtained  only  one  lime.  Any  mot  hou  that  r  aide  to 
provide  the  mode  shapes  for  displacement.  «<!«•:>»•.  h-  [>d,ng  moment .  and  shear  tor  a  biaho 
with  nonuniform  mass  and  stiffness  proper* ies  is  suitable.  I  he  approach  used  here  is  to 
calculate  the  rnodai  properties:  using  the  Tran- mission  Ma1  rix  Method  of  Martin  9 


Modal  Amplitudes 


The  modal  amplitudes  are  identified  from  blade  bending  moment  data  using  !•  .  • 
squares  fitting  approach.  The  methodology  used  here  follows  the  work  of  (laukroger  arm 
his  associates  at  the  Royal  Aircraft  Establishment  ilO,  11  referred  to  a.s  the  Xraai  Patten. 
Analysis  Method.  Their  application  of  the  method  has  been  directed  to  the  exper, mental 
determination  of  both  the  mode  shapes  and  the  modal  amplitudes,  but  the  approach  <  an 
also  be  used  when  the  mode  shapes  are  calculated  and  only  the  modal  amplitude  -  ap 
determined  from  experimental  data.  The  basis  of  the  approach  is  that  bending  mm;,- 
measurements  at  rn  stations  are  related  to  the  modal  amplitudes  of  /<  modus  b, 


{M}  -  !5|{?}  (  li) 

where  \S\  is  an  mxp  matrix.  Each  column  of  jSj  is  the  bending  moment  mode  shape  for 
the  m  stations  and  is  obtained  from  calculation  of  the  blade's  rotating  modes,  f  or  rn  ;■ 
the  modal  amplitudes  can  be  identified  by  a  least  squares  procedure 


{<?}  (;5]rl5'|)-  l\S\T{\1}  1  I 

The  matrix  product  in  Eq.  (11),  sometimes  referred  to  as  the  .S PA  matrix,  is  ;>\m  ami 
is  calculated  one  time. 


4.  APPLICATION 


The  method  proposed  here  is  examined  for  the  ca->r  of  a  <11  hi  rotor  tev.-q 
wind  tunnel  6  .  The  procedure  used  is  diagrammed  in  Figure  2.  The  spanwise  ,:ii;o  : 
measured  iri  the  wind  tunnel  are  used  to  calculate  the  blade  state  vector  ;di'-p!ace,:.<  ?  • 
slopes,  bending  moments,  and  shears)  at  200  radial  stations  using  the  mwoupled  foo-c 
response  calculation  of  |5  .  'Fhe  calculated  state  vector  is  treated  as  an  “exact"  s<  t  T 
measurements  From  this  set  of  measurements  bending  moment  “data"  are  interp, >!,.:<  d 
for  the  rn  measurement  stations.  The  modal  amplitudes  are  then  determined  !ro:n  th< 
moment  data  using  Eq.  (11)  and  the  spanwise  airloads  are  calculated  from  Eq  (s  .u  - 
(9).  The  calculated  airloads  are  then  compared  to  the  original  wind  tunnel  measurement* 
as  the  final  step.  The  upper  path  in  Figure  2  is  the  classical  forced  response  problem  w )  *  •• 
♦he  differentia!  equation  is  solved  directly.  The  lower  path  is  the  blade  hula;  <  >  •  >• 
problem  and  uses  a  modal  solution. 

fhe  CM  T1  blade  mass  and  stiffness  properties  used  for  ?h<  <  a  !<  >:  la : ,  <• , .  of  . 
modal  properties  have  been  obtained  from  12  and  are  the  same  as  tfio  e  u-  <  d  im  ‘to 
forced  response  calculation.  Twelve  flap  modes  are  c ah  ulated  at  the  nominal  rotor  d  ■  •! 
222  rpm.  It  was  necessary  to  use  31  decimal  place  accuracy  in  <  a!<  mat  mg  t  !.e  dut«  . 
of  the  Transmission  Matrix  for  the  three  highest  modus 


Modijication  to  the  Strain  Pattern  Analysis  J\feth<>d 

In  theory  it  is  possible  to  identify  the  modal  amplitudes  from  blade  displacement, 
slope,  bending  moment,  or  shear.  In  practice  it  is  difficult  to  measure  any  of  these  quanti¬ 
ties  on  a  rotor  in  flight  except  the  blade  bending  moment  and,  in  the  case  of  an  articulated 
rotor,  the  slope  at  the  blade  hinge.  The  Strain  Pattern  Analysis  Method  may  be  applied 
regardless  of  the  types  of  measurement  that  are  used.  For  instance,  if  the  root  flapping 
angle  is  to  be  used  as  a  measurement  station  then  this  is  accomplished  by  adding  an  ad¬ 
ditional  row  to  S'  using  the  root  angle  Calcuh.ted  for  each  blade  mode.  However,  if  this 
's  done,  it  is  necessary  to  normalize  (or  weight)  the  data  such  that  the  measurement  for 
'he  root  flapping  angle  is  treated  in  an  equivalent  manner  to  the  bending  moment  data. 
An  alternative  approach  is  taken  here.  The  amplitude  of  the  first  mode.  q\kr-  >s  initially 
hier.t  ihe-i  on  the  Pas  is  of  the  root  flapping  angle.  The  bending  moment  rneasu:  ements  are 
!.<••;  .  .rt :<■  i  to  account  for  the  first  mode’-,  amplitude 


t 

t 


m ,  i  m,,  •, !  n  * ,  i  :  i  'J  ] 

i.i-r*  ”  .  -  :h<  moment  mea.-uo-d  ,,t  the  tin  ,s  .  s.  is  appropriate  element  of  the 

•  :  .1  ,s  .  and  rn, ,  is  the  corr-c  ted  mea.\.nr<  motif.  The  remainder  <»f  t.;i<  mod  a: 

m.;. lit ude-  are  identified  using  hq  (11)  except  now  the  >  does  not  inr'ude  the  fir- 1  mode, 
i  ;,»•  dt.a'«  of  1 1 1 e  i'l  l  mode  amphtude  is  !;m.  corrected 


n 


(id) 


op,  u,  the  root  flapping  angle  and  nln  is  the  slope  at  the  hinge  that  has  been 
I'ci.i. for  eacn  of  the  p  rot  ating  inodes.  !  he  process  shown  ft  ere  is  iter  at  ed  upon  until 
'  ‘.e  -  fuinee  in  the  estimated  amplitude  o'  the  <::vt  mode  is  te«s  than  0. 1  •<  . 

i’  .  •  .r  - 


baseline  case  for  ‘.'is:  ^t  ’ids  av-e;*.."-  measurement  stations  and  10  modes. 
■  n  'V.ise  these  Value-,  v.u.  1  :  f,,;  <d!  i  V  rest  1 1*1  to  be  •  b  >■•  n  Figure  3 

i  -ie.r'v  1  it  h  naruionii  '  -;»a  "os  u.<  airload  itist  rif.u*  a»ti  for  t  he  lot'  k‘  <  ’.  shut* 
e  ‘  1  •  ■  d  in  f he  wind  t um.<  i  i  uuiss  dl  -t ribu; ion  for  this  blade  ■  dis.u.trmuo.:-. 


d.s*  r  :  *  'o.u  s  a  t  oust  an!  fr  tor  :  i  :  S'  and  (0).  this  u  *  s ,  ■  >  •  i  ■ 

I"  !  :  i .  ip.'  1 1  ,t  t  ♦  • 1  a;P  u<>  Iror  1  g-.r«  d  t  ! ;  a  t  :  tl»  Ulus'- 

.  c.'ong  !*.{  '  n  <  ;  1’ep.enl.  1  at  .  ••  '  1  i  r  ; !  .1  re.u-  e  : 

;  ■  t  oi  • f ■  : i].-t  r  u in*  i :*  a’  i,  ,•  .- .  .  ,  ;  I  he  air  loads  o*  <-•'>:  ,  .  pi ,  ; . • .:  .  '  a : , g, 

'  tin  ;  1 1 .  i  i  l  Ml  t  ion  fo:  "he  ret--  •  ,.'.e-  ,  U  'ht 


i  •  u 


1 1 


‘he  rei 


The  blade  mode  shapes  all  go  to  zero  at  the  blade  hinge  poin:  of  0  O.'Pdf  ]  h,.  j, 
running  mass  at  this  point  is  approximately  20  times  the  value  it  takes  on'  b«»ar<!  •..•••  .• 
the  lifting  sec tion  of  the  blade.  An  additional  consequence  of  multiply  it.g  i  In  biah--  • . - i . 
shapes  by  the  mass  distribution  is  that,  as  the  blade  hinge  approac  h«-‘.  ’?•<  no,  .  ..  • 
the  mode  shapes  and  the  mass  distribution  sometimes  result  in  a  "tad"  inboard  of  ’in 
blade  root  cutout  at  0.1GR  even  though  there  is  no  lifting  aurfe.ee  there 

The  steady  and  firs*  five  harmonies  of  airload  fee  'he  !'></  kuo* 
case  are  shown  in  Figure  -1  A  floating  scale  i.-,  used  for  :  in  -e  plots  so  ; ;  d 

the  loading  may  be  observed.  In  general,  the  agreement  between  die  <  air  if;  :  ion  u  mg  •  i 
bending  moments  and  the  measurements  is  good  and  this  is  especially  not  •<  ,t  j  , 

inboard  stations.  However,  the  method  is  unable  to  resolve  very  rapid  vana'iom  in  a. no. id 
that  are  seen  near  the  blade  tip  (see,  for  example,  the  1st  cosine  harmonic  ioad)  In  -c.tue 
cases,  for  example,  the  5th  cosine  harmonic,  disagreement  is  seen  for  nearly  all  of  the  b'ad. 
stations. 

The  airload  at  the  tip  of  the  blade  is,  of  course,  zero.  In  most  cases,  the  calculate  d 
airload  is  approaching  this  zero  airload  condition,  but  in  some  cases  (for  example  the 
3rd  sine  harmonic),  the  predicted  airload  has  converged  to  a  nonzero  value.  As  will  be 
discussed  below,  the  airload  at  the  blade  tip  is  a  measure  of  the  accuracy  of  the  sumrnai  ion 
of  Eq.  (8)  and  (9).  The  degree  to  which  a  zero  value  is  achieved  is  a  measure  of  how  vvei; 
the  fitting  of  the  measurements  is  done. 

The  agreement  between  the  measurements  and  the  calculations  art'  shown  as  a 
function  of  blade  azimuth  in  Figure  5  for  two  radial  stations,  flit'  steady  and  first  Ur 
harmonics  of  the  airload  are  shown  at  the  top  of  the  figure  while  the  vibratory  pop  ion  <>: 
the  airload  (3  to  10  harmonics)  is  shown  at  the  bottom  of  the  figure.  The  agr>  «.  meat  i-  sis  : 
to  be  quite  good  at  0.90  R ,  both  for  tin'  full  range  of  harmonic .-  and  for  <  lit'  viura'op,  par 
This  agreement  is  typical  of  all  stations  inboard  of  0.9011.  However,  a:  0.9'R  .■dgmiic • 
differences  appear  and  the  agreement  is  no  longer  satisfactory,  either  here  or  at  the  stations 
further  outboard. 

The  measurements  and  calculations  are  computed  a  ng  a  t tiree-duiteji'ioi.ai  plot 
for  the  same  150  knot,  5’  shaft  angle  cause  in  Figure  0  Tin-  method  of  pr<  .se-ntat  mo 
is  particularly  useful  in  obtaining  a  qualitative  understanding  of  *!]<•  loading  behavior  oi 
a  rotor.  The  calculated  airloads  show  very  similar  hehavim  over  mo;  t  of  t  In-  rote.:  d’>  k. 
However,  near  the  blade  tip  the  calc  ulated  loading  is  clearly  r-Oueed  compared  to  t 
measurements.  I  he  calculated  nonzero  airload  is  particularly  no* u  e.tide  wi'h  thi"  type  o! 
plot . 

Detailed  measurements  were  obtained  bv  Rabbott  ft  at.  o  f. ten  combination- 
of  air.- peed  and  sjlH|t  angle  In  general,  the  agr-ement  ' hat  hits  t«  snowi;  iot  '  r . «  1  oi 
knot,  5  shaft  angle  condition  was  observed  for  the  other  '-.si  conditions  a.s  well.  ;".i 
examples  arc  show  n  in  Figure  7  tor  t  i  <•  :■  Hoad'  at  0  91  iR  .  i  -.gup-  Tin!  -  now  t  '<•.<  1  10  k : ,  ■  > 1 
o  ,-hatt  angle  <  .ise  vhuli  •  <>rre  ponds  mo.t  1  \  to  the  f;;g\i  ••>.:  t  as.-  r  *  :>-»•*<  i  n, 

Scheir  an  3  figure  7'r.;  she -.vs  'he  75  knot.  0  shaf*  angle  «  use  v  '•  j,  app  be  * 


most  severely  loaded  case  tested  in  the  wind  tunnel.  Both  cases  show  good  agreement  for 
this  station  and  this  agreement  is  typical  of  the  inboard  stations  as  well. 

The  airload  at  the  blade  tip  must  be  zero  as  was  noted  above.  For  the  cosine 
harmonic  airload  at  the  tip 

Fzkc{  1)  --  m(l)  ^(-'n  -  k2n2)qnkc  (14) 

n  —  1 

where  all  the  mode  displacements  are  normalized  to  one  at  the  tip.  The  degree  to  which 
this  summation  approaches  zero  for  each  harmonic  airload  is  a  measure  of  how  well  the 
measurements  have  been  fitted.  The  behavior  of  tins  summation  as  the  number  of  terms 
is  increased  is  examined  in  Figure  8  for  four  of  the  harmonic  airloads  shown  previously  in 
Figure  4.  The  first  harmonic  airloads  appear  to  be  converging  to  a  nonzero  value  although 
this  is  not  completely  clear.  The  third  harmonic  airloads,  on  the  other  hand,  have  clearly 
converged  to  a  nonzero  tip  airload.  It  appears  that  only  errors  in  fitting  the  first  or 
second  modes  can  explain  the  differences  that  are  seen  here.  The  modal  amplitudes  have 
been  identified  here  using  the  root  flapping  angle  and  bending  moment  measurements. 
Although  this  provides  a  good  estimate  of  the  blade  airioads  out  to  0.90R,  the  summation 
in  Figure  8  suggests  that  if  the  known  boundary  condition  of  zero  lift  at  the  blade  tip  were 
enforc  ed  during  the  identification  process  an  improved  estimate  of  the  blade  airloads  could 
be  achieved. 

EJJe c  t  oj  N u  m  be  r  of  Modes 

Airloads  were  calculated  for  the  CK  34  rotor  assuming  20  measurement  stations 
and  using  a  variable  number  of  modes  from  4  to  12.  In  general,  the  predicted  airloads  were 
not  strongly  affected  even  when  only  four  flapping  modes  were  used.  The  3rd  harmonic 
airloads  show  the  largest  differences  for  the  150  knot,  -  5°  shaft  angle  case.  These  are 
shown  in  Figure  9  for  4,  8,  and  12  modes.  Some  improvement  is  seen  in  going  from  4  to  8 
modes,  although  even  the  four  mode  fit  gives  a  reasonable  approximation  of  the  loading. 
Little  effect  is  seen  in  going  from  8  to  12  modes  except  for  the  3rd  cosine  airload  near 
the  blade  tip  where  a  better  fit  is  achieved.  In  this  latter  case  the  modal  fit  appears  to 
show  behavior  similar  to  a  Fourier  series  representation  of  a  time  history;  that  is,  the  more 
irregular  the  behavior  the  more  terms  (modes)  that  are  required  in  the  series. 

EJJe  ft  «>/  Measurement  Error  and  Number  of  Measurement  Stations 

I  he  number  of  measurement  stations  must  be  equal  to  or  greater  than  the  number 
"f  modes  lor  the  identification  of  the  modal  amplitudes.  For  error-free  data,  there  is  little 
diffe  'reii!  »•  whether  the  number  of  measurement  stations  is  the  same  as  the  number  of  modes 
or  twice  the  number  When  errors  are  introduced  into  the  data,  however,  redundancy  in 
the  number  of  measurement  stations  shows  some  improvement  in  the  calculation  Three 
types  of  error  will  occur  in  norma!  experiment  a!  measurements  on  a  rotor  blade  of  the  kind 
di.oi  ssed  here.  The-  measurement,  at  station  r  is  expressed 

m[rj)  1  *  >  ,  ( r  j  rri(  r,  l )  *  '-(rj  ‘  f,U)  (15) 


f. 


where  rh(rj)  is  the  true  value  of  the  parameter  measured.  Th  #*  error  <  j  ( r  !  :  a  <■  ,1  i  ■ r 
or  scale  error  and  will  vary  for  each  measurement  station  along  rhe  blade.  !• 
all  harmonic  load  calculations.  The  error  ,.,(rj  represents  a  static  or  /,  ro  ob 
and  this  will  only  influence  the  calculation  of  the  steady  airload.  The  him.  <  ro,. 
represents  random  errors  distributed  in  time  and  it  is  assumed  that  these  error-  ; 
reduced  by  averaging. 

The  effect  of  calibration  errors  in  the  measurements  is  evaluated  by  assum:* 
t.5?c  errors  are  distributed  randomly  over  the  various  jnea.su rerun. ;t  >.*«!  .ons.  .-> 
shows  the  3rd  harmonic  airload  for  the  150  knot,  5  shaft  angle  rate  tor  ;*\  i 
20  measurement  stations  and  a  i  f>%  error.  Ten  modes  were  used  for  these  r.d  u; 
In  general,  the  effect  of  errors  of  this  size  was  small  and  the  3rd  cosine  harmoni*  , 
shown  here  is  typical  of  most  of  the  cases.  The  3rd  sine  harmonic  air'r-ad  sh.owt.  in 
10  represents  the  greatest  variation  that  was  seen  and  it  is  not  ibvious  that  the  <  alct 
using  the  most  measurement  stations  gives  the  best  results.  The  airload  at  0  \)0U  i-- 
in  Figure  11  for  this  same  case  and  it  is  seen  that  the  effect  of  errors  of  this  si/e  is 
These  results  show  that  the  method  used  here  is  relatively  insensitive  to  errors  th; 
bo  expected  to  occur  in  a  carefully  performed  experiment,  and  this  is  encouraging. 

Effect  of  Change s_in ..Mass  and  Stiffness  Properties 

Detailed  measurements  to  identify  the  rotating  modal  frequencies  of  a  fu 
rotor  are  riot  often  made,  but  often  the  frequencies  of  the  actual  rotor  are  lower  t  i: 
predicted  values  either  because  of  slight  mass  increases  or  stiffness  decreases  in  ecu 
tion.  or  both  To  study  •  his  effect  here  it  is  assumed  that  the  simulated  bending  r, 
data  represent  the  actual  measurements  and  the  airloads  are  calculated  using  ! 
properties  that  are  5l7  stiffer.  or  (2)  properties  with  .Vr'f  lower  mass.  These  cabm 
are  shown  in  Figure  12  for  the  150  knot,  5  shaft  angle  case  for  tin-  8th  harm-on  ■-  a 
In  most  rases,  the  **  1  for '  ■  of  mass  and  stiffness  difference--,  an  light  ns  i  seen  f.  -  -  ’ 
sm«-  airload.  In  the  ca.sc  of  the  Mti  cosine  airload,  the-  nr  ass  and  stdfness  - 
ail  of>servitbie  effect.  but  'fits  effect  is  not  large.  The  airload  <-;d<  idat'.otj  ti.e- }.  c  ;-r 
here  is.  in  general,  u  a  r.sitiv*-  to  the  small  differences  m  mas  and  si  i  fir  ess  ;,-«,(•< 
may  be  expe*  md  oetween  design  values  and  actual  Intros-..' 


5  e.M.cn.ATh  >S  <  >*■  luxi'T  SHi:.\  /i’s 


On.1  e  ' modal  amplitudes  have  been  calculated,  it  ,s  jo-  o.  t-. 

'  ti*' s t ,t ’ <•  v <■<  1 1 pre qi*  rr  e  s  il. ,r,g  t !:*■  blade,  u  it  just  :  I-'- - ' isp.ac  *  n a s  t  m •  h 
here  to  <-stiruate  the  blade  airloads  Thus,  the  same  :nethodolog\  can  he  u- 
t 1  <•  root  shear  -,  a  t.n  1.  are  important  in  understanding  th*  vir-ruf  o'*  h  ;  n e  . 

1  his  approach  or  a  umdar  one  has  be-n  used  in  the  past  In  i:  -in.  m.s  m 
identify  the  blade  root  shears  using  blad**  bending  u -•  r , ?  oat;.  I'-;  1',  \ 

t  he  modal  an. put  id*--  I  hat  have  been  us*  d  m  i  h  pa  i<  ■  '  ,  :  ,ii  l.i  •  e  1, 

can  also  f»e  n-r  <•  1  <>  es  •  irnat*-  'he  tdade  ro*>t  di;  i . .  d  ’la--  :  an  h.  .  >•; 
orig.r  a  I  --t  a  ’  •  ve<  :<■-  '  ,d  -  ie-  1  is  is  show  n  i I  ii'i.tr  i  .-,  for  t  t.e  ,  >  •  ■■  .  ■  . '  e 


\s  a- 

in,  ! 


The  effect  of  errors  in  the  data  is  examined  by  assuming  a  random  bi-'r  ; 
of  calibration  errors  at  the  various  measurement  stations.  It  is  shown  tha*  -om*  m 
men!  in  predictive  capability  is  obtained  with  measurement  redtindancv.  but.  :  g*-.-  ,  m 
the  method  is  tolerant  of  error.,  of  the  size  that  can  be  expected  in  caret, hiv  p«  norm*  •: 
experiments. 

Finally,  the  method  is  evaluated  for  differences  in  mass  and  stiffm-s-  ,.r  ■ 

might  realistically  be  expected  between  calculated  values  and  ihe  actual  rotor  prop.  :  ■ 
The  predicted  airloads  are  insensitive  to  changes  of  the  order  of  5°<  in  mass  ,,r 

The  method  proposed  here  is  shown  to  be  quite  robust  to  the  sorts  <  -f  »•>: m  r  ;  ’ 

problems  that  could  be  expected  to  occur  in  the  testing  of  a  full-scale  or  mod.  '-o  ,$;«•  :■  •  >, 
It  is  shown  that  the  computational  requirements  for  applying  the  method  arc  th"  same  <u 
similar  to  computations  that  are  normally  performed  during  th*'  design  of  a  rotor  ■  \ stem 
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Fig.  1.  Comparison  of  flap  bending  moments  for  5  ‘o  10  harmonics  on 
a  CH  34  rotor;  V  150  knots,  a*  5  .  fa)  Measured  Hap  bending 
moments  from  Reference  6.  (b)  Calc  ulated  Hap  bending  moments 

using  measured  airloads  5 


AIRLOAD,  lb/in 


Fig.  2.  Diagram  of  calculation  process. 
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•  DATA 

-  ACTUAL  MASS  DISTRIBUTION 

-  MEAN  MASS  DISTRIBUTION 


Fig.  3.  Comparison  of  calculated  and  measured  steady  (Oth  har¬ 
monic)  airload  for  CH-34  rotor;  V  =  150  knots;  a,  =  —5°. 
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Fig.  5.  Comparison  of  calculated  and  measured  time  histories  of 
airload  for  CH-34  rotor;  V  =  150  knots,  as  =  -5°.  (a)  r/R  =  0.90, 
0  to  10  harmonics,  (b)  r/R  =  0.95,  0  to  10  harmonics,  (c)  r/R  =  0.90, 
3  to  10  harmonics,  (d)  r/R  =  0.95,  3  to  10  harmonics. 
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MODE  NUMBER 

Fig.  8.  Calculated  tip  airload  with  increasing  terms  of  summation. 
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Eifc.  9.  Effect  of  the  number  of  modes  on  the  3rd  harmonic  airloads 
for  the  CH  34  rotor;  V  —  150  knots,  a,  —  •  5°. 
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Lift  10.  Effect  of  measurement  error  on  the  3rd  harmonic  airload 
for  different  numbers  of  measurement  stations  for  the  CH  34  rotor; 
V  ISO  knots,  ns  5". 
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Fig.  11.  Effect  of  measurement  error  on  time  history  at  0.90R  for 
CH-34  rotor;  V  =  150  knot,  at  =  —5°.  (a)  0  to  10  harmonics, 

(b)  3  to  10  harmonics. 


•  DATA 

-  CH-34  PROPERTIES 

-  5%  STIFFNESS  INCR. 

-  555  MASS  DECR. 


Fig.  12.  Effect  of  5%  change  in  mass  and  stiffness  properties  on  8th 
harmonic  airload  for  CH  34  rotor;  V  -  150  knots;  ft,  5J. 
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Fig.  13.  Root  flapping  shears  for  harmonics  1  to  5  calculated  from 
“measured”  bending  moment  data  compared  to  “measured”  shears 
for  CH-34  rotor. 
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lb  AdsIM'  t 

A  met  hod  is  developed  to  estimate  the  blade  normal  airloads  hv  using  him:-.'.; 
: lap  bending  moments;  that  is,  the  rotor  blade  is  used  as  a  torce  balance.  The 
blade '  rotating,  in  vacuum  modes,  are  calculated  and  the  airloads  are  then 
expre  •  >ed  as  an  algebraic  sum  of  the  mode  shapes,  modal  amp  1  i  t  tides ,  mass  ui:  -I  i  i 
imtio:.,  and  t  requenev  properties.  The  modal  amp  1  i  t  tides  are  ideiitilied  !  ror  I  in  ■ 
blade  bending  moment  s  using  the  Strain  Pattern  Analysis  Method..  The  ii'p'ic..t  i 
t  tin  met  hod  is  examined  using  simulated  flap  bending  moment  data  that  have  he 
i  i  I  cm  1  at  ed  iron  measured  airloads  for  a  till  1-scale  rotor  mi  a  wind  tunnel.  I  hi 
'"’I  im.«t  eti  lit;  a  is  are  compared  with  the  wind  tunnel  measurements.  Die  elicit 
t  he  number  •  o  me  mirement  s ,  the  number  of  modes,  md  errors  i  n  the  measurement 
Mii  t  bo  blade  proper  t  i  es  ale  examined,  and  the  method  is  shown  to  i-o  ish.i-t. 
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